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Abstract. Prospects for direct dark matter detection are 
considered in the Next-to-Minimal Supersymmetric Standard 
Model (NMSSM). The effective neutralino-quark Lagrangian is 
obtained and event rates are calculated for 73 Ge. Accelerator 
and cosmological constraints on the NMSSM are included in our 
numerical analysis, which allows us to conclude that the light- 
est neutralino from NMSSM could be detected in future dark 
matter experiments. 



1. Introduction 

The lightest supersymmetric particle (LSP), the neutralino, is considered 
now as most promising candidate, which may contribute the main fraction 
of so-called cold dark matter. 

In this paper we investigate direct detection of the relic LSP in the 
Next-To-Minimal Supersymmetric Standard Model (NMSSM) f| §. The 
NMSSM is mainly motivated by its potential to eliminate the so-called fi 
problem [§ of the Minimal Supersymmetric Standard Model (MSSM) (§, 
where the origin of the fj, parameter in the superpotential Wmssm = ^H\H-2 
is not understood. For phcnomcnological reasons it has to be of the order of 
the electroweak scale, while the "natural" mass scale would be of the order 
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of the GUT or Planck scale. This problem is evaded in the NMSSM where 
the fi term in the superpotential is dynamically generated through /i = Xx 
with a dimensionless coupling A and the vacuum expectation value x of the 
Higgs singlet. Another essential feature of the NMSSM is the fact that the 
mass bounds for the Higgs bosons and neutralinos are weakened. While in 
the MSSM experimental data imply a lower mass bound of about 20 GeV 
for the LSP §, very light or massless neutralinos and Higgs bosons are not 
excluded in the NMSSM § [§. The above arguments make an intensive 
study of the NMSSM phenomenology very desirable. Previously the Higgs 
and neutralino sectors of the NMSSM were carefully studied in []6[|— p~T[] . 
The calculation of the LSP relic abundance in the NMSSM was performed 
for the first time in Jj^] and recently in [[[3) . 

The outline of this paper is as follows. In Sec. || we describe the La- 
grangian of the NMSSM. Section || collects formulas relevant for calculation 
of the event rate for direct dark matter detection in the framework of the 
NMSSM. In Sec. | we discuss the constraints on the NMSSM parameter 
space which are used in our analysis. In Sec. |^ we describe our numerical 
procedure and discuss the results obtained. Sec. contains our conclusion. 



2. The Lagrangian of the NMSSM 

The NMSSM superpotential is g (e u = -£21 = 1): 

W = Xe i:j HiH^N ~\kN 3 
o 

+h u e VJ Q l UH 2 l - h d e tJ Q l DH( - h e s l3 L l RHi, (1) 

where Hi = (H^H~) and H 2 = {H+,H$) are the SU(2) Higgs doublets 
with hypercharge — 1/2 and 1/2 and N is the Higgs singlet with hypercharge 
0. The notation of the fermion doublets and singlets is conventional, gener- 
ation indices are omitted. The electroweak gauge-symmetry SU(2)i xU(l)y 
is spontaneously broken to the electromagnetic gauge-symmetry U(l) em 
by the Higgs VEVs (ff?) = Vi with i = 1,2 and (N) = x, where 
v = y/vf + v\ = 17 A GeV, tan/? = vijv x . 

The most general supersymmetry breaking potential has the from H 

V SQit = mllH^ + milHtf + mllNl 2 

+m 2 Q \Q\ 2 + ml\U\ 2 + m 2 D \D\ 2 + m 2 L \L\ 2 + m%\R\ 2 

~{XA x e l3 H\H 3 2 N + h.c.) - (hA k N 3 + h.c.) 

+{KA u e ij Q i UHi - h d A D e l jQ l DH{ - h^sE^U RH{ + h.c.) 

+ ^MA Q A Q + ^M'A'A'. (2) 
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The minimization conditions for the scalar potential dV/dv\^ = and 
dV/dx — eliminate three parameters of the Higgs sector which are nor- 
mally chosen to be m\, m\ and m\. 

Therefore at the weak scale we are left with the following free parame- 
ters of the NMSSM. There are the ratio of the doublet vacuum expectation 
values, tan /?, the singlet vacuum expectation value x, the couplings in the 
superpotential A and k, the parameters A\, as well as Ajj, Ajj, Ae 
(for three generations) in the supersymmetry breaking potential, the gaug- 
ino mass parameters M and M', and the scalar mass parameters for the 
squarks tuq ue, and sleptons itie,e- 

Assuming CP conservation in the Higgs sector, the Higgs CP-even and 
CP-odd 3x3 matrices are diagonalized by the real orthogonal 3x3 matrices 
II s and U p , respectively, 

Diag(m 2 5i ,m 2 52 ,miy = U sT M 2 S U S 7 
Biag(m 2 Pl ,m p2 ,0) = U pT M 2 P U P , 

where ms 1 < ms 2 < ms 3 and mp 1 < mp 2 denote the masses of the mass 
eigenstates of the neutral scalar Higgs bosons S a (a = 1,2,3) and neutral 
pseudoscalar Higgs bosons P a (a — 1,2) M. The Higgs sector of the 
NMSSM contains five physical neutral Higgs bosons, three Higgs scalars, 
two pseudoscalars, and two degenerate physical charged Higgs particles C^. 
In our numerical analysis we have included 1-loop radiative corrections to 
all Higgs mass matrices in accordance with [|l0[ |l4| . 

The masses and mixings of the neutralinos in the basis ^ T = 
(—iXi , — , ^h 2 ' * Ar ) are determined by the Lagrangian |^| 

C = -i* T A/^ + h.c. 

The mass of the neutralinos is obtained by diagonalizing the 5x5 matrix 
M with the orthogonal matrix N: 

C - — m ■ ^0 v° v° - ( X i I 

L, — — - rrit Xi Xi > Xi — I — o y 

with Xi = Nij'fyj and diag{mi,m2, ...,7715} = N M N T . The neutralinos 
Xi {i = 1-5) are ordered with increasing mass \rrii\, thus x = Xi is the LSP 
neutralino. The matrix elements M%j — 1-5) describe the composition 
of the neutralino Xi m the basis ^j. For example the bino fraction of 
the lightest neutralino is given by Af^ 1 and the singlino fraction of this 
neutralino by Af 2 5 . The remaining particle content is identical with that of 
the MSSM. 
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3. Neutralino-nucleus elastic scattering 



A dark matter event is elastic scattering of a DM neutralino from a target 
nucleus producing a nuclear recoil which can be detected by a suitable 
detector. The corresponding event rate depends on the distribution of the 
DM neutralinos in the solar vicinity and the cross section of neutralino- 
nucleus elastic scattering. 

An experimentally observable quantity is the differential event rate per 
unit mass of the target material 

dR 

&E~ r 







L m x _ 





dvf{v)v^(v,E r ), q 2 = 2M A E r . 



dq 2 

Here f(v) is the velocity distribution of neutralinos in the earth's frame 
which is usually assumed to be a Maxwellian distribution in the galactic 
frame. w max = v csc w 600 km/s and p x — 0.3 GeV-cm -3 are the escape 
velocity and the mass density of the relic neutralinos in the solar vicinity; 
fmin = (MAE r /2M 2 ed Y^ 2 with Ma and M rc( | being the mass of nucleus A 
and the reduced mass of the neutralino-nucleus system, respectively. 

The differential event rate is the most appropriate quantity for com- 
paring with the observed recoil spectrum and allows one to take properly 
into account spectral characteristics of a specific detector and to separate 
the background. However, in many cases the total event rate R integrated 
over the whole kinematical domain of the recoil energy is sufficient, ft 
is widely employed in theoretical papers |II|-[[n^|, pll pS[ ], for estimating 
the prospects for DM detection, ignoring experimental complications which 
may occur on the way. In the present paper we use the total event rate R, 
because we are going to perform a general analysis aimed at searching for 
domains with large values of the event rate R. 

A general representation of the differential cross section of neutralino- 
nucleus scattering can be given in terms of three spin-dependent J~ij(q 2 ) 
and one spin-independent Ts{q 2 ) form factors as follows |2p|| 



d(J . o\ 8Gf / 9 n~2 / 2\ t-2 / 



2\ 



a\-TU q 2 ) + c 2 -A 2 T 2 {q 2 )). 



The last term corresponding to the spin-independent scalar interaction 
gains coherent enhancement A 2 (A is the atomic weight of the nucleus 
in the reaction). The coefficients ao,i,Co do not depend on nuclear struc- 
ture but relate to the the nucleon structure and the parameters A q , C q of 
the relevant low-energy effective neutralino-quark Lagrangian Jl5| JL8| 

L cS = (A ■ Xl^l5X ■ 97^759 + jj^ -Cq-xX- Qo) , (3) 



where terms with vector and pseudo-scalar quark currents are omitted be- 
ing negligible in the case of non-relativistic DM neutralinos with typical 
velocities v w 10~ 3 c. The coefficients in Lagrangian (H) have the form: 
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^ M w P q sin20 q T 3 (Af 12 - tan6 w Nn) 
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q) 2 ™ q 2 ~ ( m x + m ?) 2 

4qL (j>qR) 



m 2 ql — (m x + m q ) 2 m? 2 — (m x + m q ) 2 



Qaii = {N12 - ts,nO w Nn)[U^ 
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<PqL = 
<f>qR = 

Pa = 



\/2AMs [U&tfu + U* 2 M 13 ] - 2v / 2/cC/ Q 5 3 AA 2 5 , 

Af 12 T 3 +Afn(Q-T 3 )tane w , 

tan 0\y Q A/11, 



n) 



A/i 
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sin/3 



■1 Tl M 3 
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The coefficients and C q take into account squark mixing — 
and the contributions of all CP-even Higgs bosons. Under the assumption 
A = k = these formulas coincide with the relevant formulas in the MSSM 
P6| . In what follows we use notations and definitions of our papers jl9) . 
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4. Constraints on the NMSSM parameter space 



Assuming that the neutralinos form a dominant part of the DM in the uni- 
verse one obtains a cosmological constraint on the neutralino relic density. 
The present lifetime of the universe is at least 10 10 years, which implies 
an upper limit on the expansion rate and correspondingly on the total 
relic abundance. The contribution of each relic particle species x nas to 
obey p2| : f2 x /ig < 1, where the relic density parameter Q x — p x j ' p c 
is the ratio of the relic neutralino mass density p x to the critical one 
p c = 1.88 • lQ- 29 hl g-cm- 3 . 

To have the papameter £l x ti^ under control we calculate it following the 
standard procedure on the basis of the approximate formula [^3| : 

«y* = 2 .i3 x io-" (|)' (^'x S? (j;^) • 

Here T 7 is the present day photon temperature, T x /T 7 is the reheating 
factor, xp = Tp/m x w 1/20, Tp is the neutralino freeze-out temperature, 
and Np is the total number of degrees of freedom at Tp. The coefficients 
a,b are determined from the non-relativistic expansion (a ann v) ~ a + bx 
of the thermally averaged cross section of neutralino annihilation in the 
NMSSM. We adopt an approximate treatment not taking into account 
complications, which occur when the expansion fails [p4| . We take into 
account all possible channels of the \~X annihilation in the NMSSM Ji3| ]. 

Since the neutralinos are mixtures of gauginos, higgsinos, and singlino 
the annihilation can occur both, via s-channel exchange of the Z a and Higgs 
bosons and t-channel exchange of a scalar particle, like a selectron. This 
constrains the parameter space, as discussed by many groups [^3[ p5| . 

In the analysis we ignore possible rescaling of the local neutralino den- 
sity p which may occur in the region of the NMSSM parameter space where 
fi x /ig < 0.025 p6[-p^]. If the neutralino is accepted as a dominant part of 
the DM its density has to exceed the quoted limiting value 0.025. Otherwise 
the presence of additional DM components should be taken into account 
and the SUSY solution of the DM problem with such low neutralino density 
becomes questionable. We assume neutralinos to be a dominant compo- 
nent of the DM halo of our galaxy with a density p x — 0.3 GeV-cm -3 in 
the solar vicinity and disregard in the analysis points with f2 x /ig < 0.025. 

The parameter space of the NMSSM is constrained by the non- 
observation of the supersymmetry particles with high energy colliders LEP 
and Tevatron. Due to the possible presence of singlino component, the 
Zx°X° an d ZS a P a couplings can be suppressed in the NMSSM. Therefore 
the neutralino and Higgs mass bounds are much weaker than in the MSSM 
[||. The consequences from the negative neutralino search at LEP for the 
parameter space and the neutralino masses have been studied in |g|, ||. 
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We used the following constraints from LEP. For new physics contribut- 
ing to the total Z width Ar(Z -> x+x~ + Z -> x°X°) < 23 MeV. For new 
physics contributing to the invisible Z width AT(Z — » X^Xj) < 8 MeV. 
From the direct neutralino search B(Z — > x?Xj) < 2 x 10~ 5 for j = 2, . . . , 5, 
and B(Z — > x° X°) < 5 x 10~ 5 , for i,j = 2, . . . , 5. We have also included the 
experimental bounds from the direct search for pseudoscalar Higgs bosons 
produced together with a Higgs scalar at LEP , but this does not signif- 
icantly affect the excluded parameter domain S. 

In our numerical analysis we use the following experimental restric- 
tions for the SUSY particle spectrum in the NMSSM: m~+ > 90 GeV for 
charginos, > 80 GeV for sneutrinos, mg fl > 80 GeV for selectrons, 
rriq > 150 GeV for squarks, > 60 GeV for light stop, m H + > 65 GeV 
for charged Higgses and ms 1 > 1 GeV for the light scalar neutral Higgs. 
In fact, it appeared that all above-mentioned constraints do not allow ms 1 
to be smaller then 20 GeV. 



5. Numerical Analysis 

Randomly scanned parameters of the NMSSM at the Fermi scale are the 
following: the gaugino mass parameters M' and M, the ratio of the doublet 
vacuum expectation values, tan /3, the singlet vacuum expectation value x, 
the couplings in the superpotential A and k, squared squark mass param- 
eters tUq i 2 for the first two generations and toq^ for the third one, the 
parameters A\ , Ak, as well as A t for the third generation. The parameters 
are varied in the intervals given below 
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For simplicity the other sfermion mass parameters m^ i2 , rn 2 Di2 , m L , 
m 2 Ei 2 , and mf^, m 2 B , Tn\ a , rn 2 E3 , are chosen to be equal to and m^, 
respectively. Therefore masses of the sfermions in the same generation differ 
only due to the D-term contribution. Other parameters (except A t ) of the 
supersymmetry breaking potential Ajj, A E , A E (for all three generations) 
are fixed to be zero. 
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Figure 1. Total event rate R for 73 Ge, the LSP gaugino fraction (Affi + 
N% 2 ) an d singlino fraction (J\f? 5 ) versus the LSP mass (from up to down). 
The left (right) panel presents results obtained without (with) taking into 
account the cosmological relic density constraint. 
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The main results of our scan are presented in Figs. 1-4 in the form of 
scatter plots. Given in Fig. 1 are the total event rate R for 73 Ge, and the 
LSP gaugino fraction (Af^ + N% 2 ) an d singlino fraction (N% 5 ) versus the 
LSP mass. The left panel in Fig. 1 presents the above-mentioned observ- 
ables obtained without taking into account the cosmological relic density 
constraint. In this case the total expected event rate R reaches values up 
to about 50 events per day and per 1 kg of the 73 Ge isotope. As on can 
see from Fig. 1 the small-mass LSP (less then about 100 GeV) are mostly 
gauginos, with very small admixture of the singlino component. Large 
masses of the LSP (larger then 100 GeV) correspond to sizable gaugino 
and singlino components together with some higgsino fraction. The results 
of implementation of the cosmological constraint 

0.025 < Q x hl < 1. 

one can see in the right panel of Fig. 1. There is approximately a 5-fold 
reduction of the number of the points which fulfill all restrictions in this 
case. Nevertheless quite large values of event rate R (above 1 event/day/kg) 
still survive the cosmological constraint. The gaugino component becomes 
more significant, but the singlino fraction can not be completely ruled out 
especially for large masses of the LSP. The higgsino component of the LSP 
remains still possible only for LSP masses in the vicinity of M z . The lower 
bound for the mass of the LSP now becomes about 5-7 GeV (Fig. 2). 
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Figure 2. Total event rate R for 73 Ge versus the LSP mass in the region 
of low LSP mass. 



Figure 3 shows correlations between mass of the LSP and mass of the 
lightest scalar Higgs boson Si. The upper (about 160 GeV) and lower 
(about 20 GeV) bound for the mass of Si are clear seen. Given in Fig. 4 
are the total event rate R for 73 Ge versus the masses of all scalar Higgs 
bosons 51,2,3- 
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Figure 3. Correlations between mass of the LSP and mass of the lightest 
scalar Higgs boson Si. 

6. Conclusion 



In not too far future new very sensitive dark matter (DM) detectors 




[ pl[ may start to operate, and one expects new, very important data from 
these experiments. The prospects of the direct and indirect detections of 
the LSP have comprehensively been investigated |jl5| in the various versions 
of the MSSM. 

In the paper we address the question whether the Next-to-Minimal Su- 
persymmetric Standard Model can be attractive from the point of view 
of the direct detection of the neutralinos provided the neutralino is the 
stable LSP. To answer the question we derived the effective low energy 
neutralino-quark Lagrangian, which takes into account the contributions of 
extra scalar Higgs boson and extra neutralino. On this basis we calculated 
the total direct-dark-matter-detection event rate in 73 Ge as a representa- 
tive isotope which is interesting for construction of a realistic dark matter 
detector. We analyzed the NMSSM taking into account the available ac- 
celerator and cosmological constraints by means of random scan of the 
NMSSM parameter space at the Fermi scale. We demonstrated that the 
cosmological constraint does not rule out domains in the parameter space 
which correspond to quite sizable event rate in a germanium detector. 
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Figure 4. Total event rate R for 73 Ge versus versus the masses of scalar 
Higgs bosons 51,2,3 (from up to down). 



Due to relaxation of the gaugino unification condition wc found domains 
in the parameter space where lightest neutralinos have quite small masses 
(about 7 GeV), acceptable relic abundance and sufficiently large expected 
event rate for direct detection with a 73 Ge-detector. 

Therefore the NMSSM looks not worse then the MSSM from the point 
of view of direct dark matter detection. The question arises: Is it possible to 
distinguish MSSM and NMSSM by means of direct dark matter detection 
of LSP? It is a problem to be solved in future. The question can disappear 
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by itself if negative search for light Higgs with LHC rules out the MSSM. 
In fact, the upper tree level mass bound for the lightest Higgs scalar of the 
MSSM 

ml < m% cos 2 2/3 (4) 

is increased in the NMSSM to m% 1 < m| cos 2 2/3 + X 2 (vf + «§) sin 2 2(3 and 
the NMSSM may still remain a viable model when the MSSM can be ruled 
out due to (Q). Therefore the NMSSM might remain a viable theoretical 
background for direct dark matter search for relic neutralinos in the post- 
MSSM time. 
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